Optical manifestations of symmetry breaking in bilayer graphene 
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We propose a spectroscopic method of identifying broken symmetry states of bilayer graphene. 
We demonstrate theoretically that, in contrast to gapped states, a strained bilayer crystal or nematic 
phase of the electronic liquid are distinguishable by the dependence of the lineshape of absorption 
on the polarization of the light. This property is characteristic for both the infrared and far-infrared 
spectral ranges, which correspond to the absorption by transitions between low-energy bands and 
split bands, and transitions between the low-energy valence and conduction bands, respectively. 



In this Rapid Communication we study how symme- 
try breaking in bilayer graphene (BLG)A manifests it- 
self in the optical spectra in the infrared (IR) and far- 
infrared (FIR) spectral ranges^ - — Symmetry breaking 
in BLG may be caused both by external perturbations 
and by internally developed instabilities generated by 
the electron-electron interactions. For example, strain 
in BLG which might be inflicted on the crystal involun- 
tarily upon thermal annealing and cooling of suspended 
BLG devices would asymmetrically change the topol- 
ogy of the low-energy dispersion^ - — Also, by apply- 
ing a perpendicular electric field, one breaks the inver- 
sion symmetry of the lattice opening an externally tun- 
able gap between the valence and conduction bandsi Al- 
ternatively, at low temperatures, undoped pristine BLG 
may undergo a spontaneous symmetry breaking transi- 
tion into one of the recently discussed strongly corre- 
lated ground statesJ^r— In particular, the phases fa- 
vored by the renormalization of short-range electron- 
electron interaction constants 2 ^ - — are a nematic state 
in which the isotropy of the band structure is reduced in 
a similar way to strained BLG ) 10 i 25 and gapped layer- 
antiferromagnetio^S, and spin flux phases^ Although 
several transport experiments 2 ^ - — reported observations 
of some broken symmetry states in BLG, the exact na- 
ture of the ground state still remains to be established. 

Here, we show how infrared (IR) and far-infrared (FIR) 
absorption spectroscopy can be used to distinguish be- 
tween some of the broken symmetries. The feature dis- 
cussed below is that strain (or a phase transition to the 
nematic state) induces a dependence of absorption of 
light on the polarization of the radiation. This anisotropy 
can be characterized by a factor 
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where g« (g±) is the absorption coefficient of light with 
linear polarization e parallel (perpendicular) to the prin- 
cipal strain axis (or the direction chosen by the order pa- 
rameter of the nematic phase). In contrast, the gapped 
phases show isotropic absorption and notable qualitative 
differences in the lineshape of their absorption spectrum 



as compared to the strained (nematic) and unperturbed 
BLG states. 

The absorption coefficient 2 ^ analyzed in this study, 
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is given by the ratio of the Joule heating and the flux 
of incident radiation. Equation ([2]) describes transitions 
between initial and final plane wave states marked by 
indices A which include the band, branch, and valley (spin 
is also taken into account); e p \ is the energy of a plane 
wave with momentum p, / is the Fermi function, A is 
the normalization area, and 



M$ = (pA|jl|pA')( P A'fo|pA) 



(3) 



are determined by the current density operator j = ed p rl. 
Since the momentum transferred by light is negligibly 
small, in the momenta of electrons in the initial 

and final states of the inter-band transitions are taken 
to be equal. The corresponding plane wave state wave 
functions are the four-component eigenstates |...) oc 
(ipAi ,£,4>B2, 4'A2 1 £,ipBi ) of the HamiltonianiiiS 
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where -Ai(2) and -Bi(2) identify the sublattices of the hon- 
eycomb lattices in the top (bottom) layers, £ — ± dis- 
tinguishes the K and K' corners of the hexagonal Bril- 
louin zone, and Pauli matrices a x , cr y , a z , and r z act 
on the sublattice and valley components of | . . .), respec- 
tively. Also, 71 stands for the inter-layer coupling be- 
tween atoms on A2 and B\ lattice sites, and v is the 
Dirac velocity for monolayer graphene. 

For unperturbed BLG (SH = 0), Eq. (g]) determines^ a 
pair of low-energy bands near the Brillouin zone corners 
with spectrum e p ± — ±p 2 /2m* = ±e , m* = 7i/2t> 2 
0.035m e for states located mainly on sublattices A\ and 
£?2- Equation ^ also determines two split bands with 
quadratic dispersion e p ± ±(71+60) and wave functions 
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that have equal weight on lattices A 2 and B\ . Note that 
the split band states are almost unperturbed by the strain 
or the formation of a nematic phase or one of the gapped 
phases X 

For the sake of convenience, the "valley momentum" 
p in Eqs. <j3j) and (j4} is determined with respect to the 
position of the Dirac point in the graphene monolayer, 
which is shifted from the Brillouin zone corners K and 
K' in a homogeneously strained crystal^ In monolayer 
graphene, such a shift p — > p' + a = p in the momentum 
space is trivially absorbed into a gauge transformation so 
that it does not influence observable characteristics, such 
as the absorption spectrum. In contrast, in BLG the in- 
terplay of such a shift and the interlayer skew hopping 
73 which couples sublattices A\ and B2 generates a per- 
turbation which cannot be eliminated from the Hamil- 
tonian by any gauge transformation. This, as well as 
other possible symmetry-breaking perturbations in BLG 
arc included in 
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The first term in SH accounts for the effect of strain^ 
or nematic order, with the direction of the unit vector 
A = (cos 29, sin 29) set by the direction £ — (cos 9, sin 9) 
of the principal stretching direction. This perturbation 
changes the low-energy electron dispersion into 
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which features a Lifshitz transition at e = ±|iu|, from 
an almost parabolic dispersion at |e| 3> \w\ to a pair of 
Dirac cones at p = ±\ / 2mw£ shifted from each other 
along the anisotropy axis £. In the following, the axis 
£ will be used as the reference direction to distinguish 
between "parallel", e || £, and "perpendicular", e _L £, 
polarizations of light and the corresponding absorption 
coefficients g\\ and g± in Eq. (TTJ). 

Also, in SH, antiferromagnetic (AF) and spin-flux (SF) 
states are characterized by a splitting 2u between the va- 
lence and conduction bands at the K point and have iden- 
tical spectra e p ± = ±\/u 2 + Cq, the same as in graphene 
with an asymmetry gap opened, for example, by a per- 
pendicular electric field.— For magnetic phases of BLG, 
S characterizes the spin quantization axis, but since the 
current operator j does not depend on the spin at all, 
these BLG states will behave identically in the optical 
absorption. 

In the IR spectral range, cj « 71 « 0.4eV, where u is 
the energy of the incident light, the optical transitions 
which are sensitive to the BLG symmetry breaking are 
those between the small momentum parts of the split 
bands and the low-energy bands. In addition, the spec- 
tral density of IR absorption includes the contribution 
from the transitions between the two low-energy bands 
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FIG. 1. (a) Absorption coefficients g\\ and g± and absorption 
anisotropy Q for undoped strained/nematic bilayer graphene 
with w — 5meV and T = 3K. (b) Sketch illustrating optical 
transitions from the low-energy bands at g, = and T — in 
strained BLG or the nematic phase of BLG. Solid (dashed) 
lines denote filled (empty) states. The two extremal parts of 
the dispersion are shown for p |] L (with Dirac points at ±p w 
with p w — \/2mw), and p _L £. The arrows mark the range of 
p for which the threshold transitions with ui — 71 + w exist, 
(c) The absorption coefficient for undoped, gapped bilayer 
graphene. (d) The absorption coefficient for the unperturbed 
system with finite chemical potential g. (e) Absorption coef- 
ficient and (f) anisotropy Q of strained BLG with finite dop- 
ing for w — 5meV and T = 3K, for g — 2meV < w and 
g — 6meV > w. 
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at high momentum which provides an almost-constant 
background in the absorption spectrum^i 3 - In Fig. [TJa) 
we show the calculated absorption spectrum for the two 
characteristic polarizations of light and the polarization 
factor Q for undoped, strained BLG with w = 5meV at 
T = 3K. The absorption spectrum of BLG in the nematic 
phase would have the same features: the anisotropy of 
the absorption for to ss 71 + w, where the absorption for 
e _L £ (gray line) shows a strong peak near the threshold 
whereas absorption for e || £ (black line) is weak: 
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This occurs because the matrix elements M in Eq. (|3J) 
preferentially select transitions with p _L £ for g\\ and p | 
£ for gj_. In the first case, when p _L £, only transitions 
from right at the center of the K point have the threshold 
energy 71+iu, but in contrast, when p || £, there is a finite 
range of momenta, — \/2mw < p < \/2mw, where the 
low-energy valence band and conduction split band are 
shifted on the energy scale by exactly 71 -I- w, [Fig. Ifb)] . 
(This is also true for the valence split band and the low- 
energy conduction band). This produces a singularity at 
the interband absorption edge, 71 + w£L 

For comparison, in Fig. HJc) we show the absorption 
spectrum characteristic for any of the gapped states of 
BLG. Here, the absorption coefficient does not depend 
on the polarization of the photon and there is a fea- 
ture at the threshold ui = 71 + u of the lowest energy 
inter-band transitions. The height of this peak is con- 
stant for u ^> ksT. When fc^T > u (this situation 
is considered having in mind the gapped state caused 
by the inter-layer asymmetry due to external perturba- 
tion rather than an intrinsic phase transition) thermal 
occupation of the low-energy conduction band also allows 
transitions to the split band with energy 71 — u, which 
yields a small additional polarization-independent peak 
in g(ui). The absorption of IR light by unperturbed BLG 
with finite doping and chemical potential /1 is shown in 
Fig. [Ud), 3 — in precise agreement with previous calcula- 
tions of the optical conductivity^' 3 -' 3 ^ 

Figure [He) shows the absorption spectrum of strained 
BLG with finite doping (the chemical potential /i 7^ is 
counted from the Dirac point energy) . Note that the ne- 
matic phase is not expected to survive at finite doping. 
When \fi\ < w (dashed line), the absorption spectrum 
remains almost unchanged as compared to the undoped 
case, but for \fx\ > w (solid line) a new peak appears 
at lo = 71 — w in the e || £ polarization but not in 
the e _L £ polarization. This occurs because, with this 
level of doping, transitions from the low-energy conduc- 
tion band above the Lifshitz transition become accessible. 
This asymmetry manifests itself in the polarization de- 
gree Q of the absorption, shown in Fig. [IJf) for \fj,\ < w 
(dashed line) and |/z| > w (solid line). 

We now turn our attention to absorption in the FIR 
frequency range, \u)\ ~ 2w <C 71, where the relevant 
optical transitions occur between the two low-energy 
bands. Figure HJa) illustrates features of the absorption 
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FIG. 2. Absorption coefficient in the FIR frequency range 
for (a) undoped strained BLG (or nematic state) with w = 
5meV, T = 3K; (b) the doped, strained case with w = 5meV, 
T = 3K, and [i = 2meV and 6meV. (c) The undoped gapped 
phase, (d) The unperturbed case with n > 0. 



by strained BLG for the two characteristic polarizations 
of FIR radiation, e || £ and e _L £: a weak polarization 
dependence described by Q ps +0.3 at energies lu « 2w 
indicating that the absorption is strongest for light po- 
larized in the direction of the principal strain axis. Note 
that for FIR light the relation between absorption in dif- 
ferent polarizations is opposite to what we found for the 
IR spectral range. A dip in the absorption at very low 
energies is due to the finite temperature. The absorption 
by strained BLG with finite doping is shown in Fig.^b). 
The only effect of the doping is to cut off the absorp- 
tion for u> < The absorption by the gapped phase 
is shown in Fig. [He) for comparison: it has no polariza- 
tion dependence but has a peak at uj = 2u corresponding 
to the threshold of the optical transition. The doped 
but unperturbed BLG, Fig.^d) shows a step in g(ui) at 
u) = 2/i£ Therefore, one can distinguish the type of sym- 
metry breaking in BLG using FIR spectroscopy through 
the weak absorption anisotropy of the strained/ncmatic 
state and the band-edge peak in the gapped state. 

In conclusion, we have shown that IR and FIR absorp- 
tion spectroscopy can distinguish between the gapless 
anisotropic nematic (or strain-induced) and the isotropic 
gapped broken symmetry states of BLG. The former has 
a characteristic strong dependence on the orientation of 
the polarization of the incident radiation, such that in 
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the IR frequency range, light polarized perpendicularly 
to the strain axis (or symmetry-breaking axis in the ne- 
matic phase) will be absorbed very strongly at u> « ji+w 
in the undoped system, whereas light polarized parallel 
to this axis acquires a characteristic feature in absorp- 
tion when the doping is such that fi > w. There is also 
a weak polarization dependence for strained BLG in the 
FIR regime, with the parallel polarization being absorbed 
more strongly. In contrast, the isotropic gapped phases 



show no absorption anisotropy, but do have qualitatively 
different lineshape from both strained BLG and the un- 
perturbed state at zero or finite doping. 
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For strain^ w = j3(5—6') 
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where S and 5' 



V dR AB 

are eigenvalues of the strain tensor, and Rab is the closest 
neighbor distance between carbon atoms in the monolayer. 
Since this divergence is created by the existence of parallel 
bands, it may be cut off by any additional perturbation 
which modifies the slope of one band with respect to the 
other. For example, the next-nearest neighbor inter-layer 
hops, the inter-layer hops from dimerized to un-dimerized 
lattice sites, or a momentum dependence of 71 (and there- 
fore of the effective mass) can all cause this effect. 
For doped BLG, there is a step at uj — 71 ± 2[i because, 
for 71 < oj < 2/j,, the transitions between the split valence 
band and low-energy conduction band are Pauli-blocked. 
There is also a very sharp spike caused by the transitions 
between the two parallel conduction bands (low-energy 
band and split band at 71) for k < ^/2m/i/fi. 
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